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1 Introduction 



Ince studied, in his well-known book of differential equations [1], the following equation 

w" + 3ww' + w 3 = q(z) 

and proved that it has the general solution w = u' /u, where u is a general solution of the linear 
equation of the third order u'" = q{z)u. This equation was also studied by Davis in [2] as a particular 
case of the generalized Riccati equations (according to Davis the family of these nonlinear equations 
was first studied by E. Vessiot in 1895 and G. Vallenberg in 1899; see [3, 4, 5] for some more recent 
studies related with higher-order Riccati equations). Later on Leach et al [6, 7] consider the equation 

q + qq + (3q 3 = (1) 

and point out that "for [3 = 1/9 is linearizable, possesses eight symmetries and is completely inte- 
grable" and they add "consequently, we could expect that this remarkable mathematical property 
corresponds to an important physical one appearing (or disappearing) for this value which conse- 
quently would appear as a critical one". This particular (5 = 1/9 equation was also obtained in [8] 
in the study of nonlinear equations with the maximum number of symmetries (see [9], [10], [11], 
and [12] for the Lie symmetry approach to dynamical systems). 

Recently Chandrasekar et al [13] have studied a generalization of this equation obtained as a 
particular case of the Lienard equation 

x + f(x)x + g(x) = 

given by f(x) = kx and g(x) = (l/9)kx 3 + Xx. Although this new equation also belongs to the 
generalized Riccati family studied by Davis and Leach et al, they make use of a two step procedure 
to solve the problem: firstly they use the so-called Prelle-Singer method [14, 15, 16, 17] for obtaining 
a set of time-dependent integrals of motion and secondly they use this time-dependent family in 
order to compute the solution. The result is interpreted, when A > 0, as an "unusual Lienard type 
oscillator with properties of a linear harmonic oscillator" . But we wish to call the attention to one 
property discussed in the final part of the paper (after finalizing with the Prelle-Singer method): 
the existence of a Lagrangian description. 

The main objective of this article is to develop a deeper analysis of these nonlinear equations using 
the Lagrangian formalism as an approach. In fact the starting point of our approach is the fact that 
the Riccati equation belongs to a family of nonlinear equations admitting a Lagrangian description. 
This has interesting consequences, the most important of them is that Riccati systems are systems 
endowed with a preserved energy function. We study the two nonlinear systems firstly in one 
dimension and then in two dimensions. Moreover we prove that the two-dimensional extensions 
are not only integrable but also super-integrable. We note that this situation has certain similarity 
with the one-dimensional nonlinear oscillator studied by Mathews et al in [18] that has been proved 
to admit a superintegrable two-dimensional version [19]. 



2 



The plan of the article is as follows: In Sec. 2 we present a Lagrangian approach to a family 
of nonlinear equations that includes the second-order Riccati equation as a particular case. Sec. 
3, that is devoted to the first nonlinear system ('dissipative'-looking system), is divided in three 
parts corresponding to the one-dimensional system, geometric formalism and symmetries, and two- 
dimensional system and super-integrability, respectively. Sec. 4, that is devoted to the second 
nonlinear system ('nonlinear oscillator'), also firstly studies the n = 1 system and then the two- 
dimensional system that, as we have pointed out, is also endowed with super-integrability. Finally 
in Sec. 5 we make some comments. 



2 Lagrangian formalism and second-order Riccati equa- 
tions 

In this article we shall consider the following nonlinear second-order equation 

y" + [6o(t) + b!(t) y y + ao(t) + ai (t)y + a 2 (t)y 2 + a 3 {t)y 3 = , (2) 
where we suppose that 03 > and the two functions bo, bi, are not independent but satisfy 

V«3 2a 3 



The more important property of this equation is that it can be transformed into a third-order linear 
equation by the substitution 

1 v'(t) 

Thus the equation (2), that is the natural second-order generalization of the well known Riccati 
equation, is therefore a nonlinear equation the solution of which can be expressed in terms of 
solutions of a linear equation of third order. In fact it can be considered as the particular n = 2 
case of a more general situation that can be approached by Lie theory or by the action of an 
operator R. From the Lie theory of symmetries of differential equations, the invariance of the 
n-order linear equation, 

u (B) +pi(t)v^ + ...+p n {t)v = 0, 

under the vector field X = v (d/dv), that represents the infinitesimal generator of dilations, means 
the existence of a change v = e u such that X becomes X = d/du, and the transformed equation 
reduces to an equation of order (n — 1) for y = u' that for the n = 2 reduces to the usual first-order 
Riccati equation. Alternatively if R denotes the following differential operator 

R = j t +y(t), 

then the Riccati equation of order n is given by 

(R n + Pl R n - 1 + ...+ Pn+1 R + Pn )y + Pn+1 = 
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where pj = Pj(t), j = 1,2, ...,n + 1, are n + 1 arbitrary functions. We note that these two 
methods lead to the subfamily of the nonlinear equations with the coefficient of the higher power 
equal to one; nevertheless the general Riccati equation is generated by the change t = /(r) of the 
independent variable; for example for n = 2 we obtain 

y" + [pi + 3y] y ' + p 3 + P2V + Piy 2 + y 3 = o 

and after the time reparametrization we arrive to 

d 2 y , f\ , ow i^ 



dr 2 



if pi + afv ± + (/' 2 p 3 ) + (7' 2 P 2 ) y + (/'V) y 2 + /' 2 y 3 = o . 



dr 



We are interested in the study of nonlinear systems given by (2) but first, in this section, we 
consider a more general familly from which the second-order Riccati equation appears as a particular 
case. 

At this point we recall that a Lagrangian function L is called 'natural' or 'of mechanical type' 
when it is of the form L = T — V, where T is a quadratic kinetic term and V is a potential 
function. Most of the known Lagrangian equations arise from Lagrangians of this particular type; 
nevertheless the Lagrangian formalism is well defined, not only for these specific functions but also 
for more general Lagrangian functions. 

Proposition 1 The nonlinear second-order Riccati equation admits a Lagrangian description. 
Proof: We first consider the following one degree of freedom Lagrangian 

v x + kU(x,t) ' V ; 



Then we arrive to the following second-order nonlinear equation 

^ + (§)*!£(§) + {^UU' X + kill = . (4) 
In the particular case of the function U = U(x,t) being a quadratic function 

U = c (t) + Cl (t)x + c 2 {t)x 2 , 

the above equation (4) reduces to 

d 2 x ., , ,/dx\ 2 o 

+ {oq + bix)y—j + ao + a\x + a 2 x + a%x =0 (5) 

where the four functions ao, ai, 02, as, are given by 

a = (^)c ci + c , ai = c c 2 + (^)c? + ci , 
3 

02 = (2)°!^ + c 2 , «3 = cl , 



4 



and the two functions 60 and 61 satisfy the appropriate restrictions 

3 

b o = {-j)ci , h = 3c 2 . 

Thus the second-order Riccati equation (5), that is a particular case of the equation (4), is the Euler- 
Lagrange equation of the Lagrangian function (3) in the particular case of a quadratic function 
U. 

As a corollary of this proposition we can state that, if the function U is time-independent, then 
the nonlinear equation (4) has a first-integral that can be interpreted as a preserved energy. The 
idea is as follows: if we restrict our study to the case of time-independent systems, that is, to 
nonlinear equations arising from a Lagrangian of the form 

(6) 



v x + k U(x) 

then we can define an associated Lagrangian energy -Ex by the usual procedure 

E L = A(L) — L , A = v x , 

dv x 

and we arrive to 

_ - (2v x +kU(x)) d 
L {v x + kU(x)) 2 ' dt L 

Note that L is non-natural and, as there is neither kinetic term T nor potential function V, the 
energy can not be of the form El = T + V. Note also that a Lagrangian is defined up to certain 
ambiguities; that is, L = c\L + cq determines the same differential equation but it leads to an 
slightly different energy El given by El = c\El — cq. In the 'natural' case, c\ is determined by the 
corresponding Riemannian metric (hence the classical one half coefficient) and cq is absorbed in 
the potential; here we have just taken c\ = 1, c$ = 0. Concerning the negative sign, that could be 
considered as something inconvenient, it does not matter at all (it is only an aesthetic question); 
in fact it can be removed just by chosen L = —L as a new Lagrangian. 

As for the 'natural' case we can obtain the solution of the dynamics from the conservation law 
of the energy. If we assume that El takes the constant value El = E, then we arrive to 

E v 2 x + 2(1 + jfe E U(x))v x + k(l + kE U{x))U{x) = 

and on solving for v x and making separation of variables we arrive at 

r Edx 

° ~ ~ lx (1 + kE U(x)) ± V 1 + k E U{x) ' 

The motion is confined to the region where E> — \/{kU). To sum up, a time-independent system 
describedby the nonlinear equation (4) is solvable and the solution of the dynamics is given (up to 
one integration) by (7). 
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An remarkable property is that the conservation of the energy El leads to a plus/minus sign in 
the expression for the velocity 

dx _ -(l + kEU{x)) ± y/l + kEU(x) 
~dt ~ E ' 

Thus, we obtain two different values for the velocity at the same point x. It is known that in the 
standard case of a particle in a potential V(x) we also have two possibilities but both with same 
modulus (the positive value for the motion from left to right and the negative for the opposite 
motion from right to left). In this case the situation is different; we have a Lagrangian but not a 
potential and the two possible values differ, not just in the sign, but in the absolute value; thus the 
motions from left to right and from right to left take place at different velocities. 

We close this section with the problem of the existence of alternative Lagrangians. 

In differential geometric terms a time-independent Lagrangian function L determines an exact 
two- form ujl defined as 

Ol = (q^-) dx ' u L = - d6 L , 

in such a way that, when L is nonsingular, ujl is symplectic and the dynamics is given by the vector 
field solution of the equation 

i(T L )uj L = dE L . 

In this particular case u>l and Tl are given by 

2dxf\dv x ^ d „ d „ , 1. , . ,„, .. rr , 

UL= (v x + kU( X ))* > rL = V -te +F -d£> F* = -(- 2 )H3v x + kU(x))U x . 

An important property of the Lagrangian formalism is that for one degree of freedom systems 
there exist many different equivalent Lagrangians [20, 21]. A sketch of the proof is as follows: in a 
two-dimensional manifold all the symplectic forms must be proportional. Hence, for a one degree of 
freedom Lagrangian, any other symplectic form oj 2 must be proportional to ul, that is uj 2 = fu>L- 
Then 

i(T L ) uj 2 = f i(T L ) uj l = f dE L 

The right-hand side is an exact one-form if, and only if, dfAdE^ = 0, which shows that / must be 
a function of El. In this case it can be proved that the new symplectic form uj 2 is derivable from 
an alternative Lagragian L 2 ^L for Tl- 

In this particular case, starting with the Lagrangian (6) and assuming for the constant of motion 
/ the particular expression / = (— 1/Ex) 3 / 2 , we have obtained the following function 

L 2 = ^2v x + kU{x) (8) 

as a new alternative Lagrangian for the t-independent version of the equation (4). This new 
Lagrangian, that is neither 'natural' or 'of mechanical type', is equivalent to L in the sense that 
both determine the same dynamics. Nevertheless we must say that it is not clear whether L 2 will 
lead to simpler expressions for other dynamical properties; so, in the following, we only use the 
original Lagrangian (6). 
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3 Lagrangian conservative approach to a 'dissipative'- 
looking nonlinear system 



3.1 One-dimensional nonlinear system: energy and integrability 

We now apply the formalism introduced in Sec. 2 to the study of the following nonlinear equation 



It is a special case of equation (2) and because of this is a Lagrangian equation. In fact, it is easy 
to verify that it can be obtained from the following Lagrangian function 



Two important properties are: 

(i) Firstly we note that the value (3 = 1/9, pointed out by Leach et al [7] as the particular value 
introducing a high degree of regularity in the nonlinear problem (1), appears now as related with 
the Lagrangian origin of the equation. That is, only if (3 = 1/9 the equation (1) belongs to the 
Lagrangian family (9) arising from (10). 

(ii) Secondly this equation looks like a dissipative equation with the damping term proportional to 
xv x ; nevertheless it is in fact a conservative system because of its Lagrangian origin. What happens 
is that the term conservative is usually considered in the Newtonian sense, that is, a particle moving 
in a one-dimensional potential and forces determined as the gradient of the potential. Here the 
force is a velocity-dependent force and conservative just means the existence of a preserved (non- 
Newtonian) energy function that is given by 



Next we turn to the solution of the dynamics. Instead of considering directly the nonlinear 
equation we can solve this problem by making use of the conservation law of the energy; if we 
assume El = E, then we arrive at 




(9) 



E L = 



- (2v x + kx 2 ) 
(v x + k x 2 ) 2 



Evl + 2(l + kEx 2 )v x + k{l + kEx 2 )x 2 = 



and solving for v x we obtain 



dx 
~dt 



(1 + k Ex 2 ) ± V 1 + k Ex 2 
E 



So, after integration, we arrive at 



1 

kx 



(1 ± v 7 ! + k Ex 2 ) 
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which yields 

2t 

X ~ kt 2 -E 

that represents the solution of the dynamics as a function of the constant value E of the Lagrangian 
energy (for easy of notation we give the solution for the particular initial conditions (to = 0, xq = 0)). 
The system is well defined for (k > 0, E < 0) and (k < 0,E > 0), but for (k < 0, E < 0) or 
(k > 0, E > 0) it is singular at t = ±yf\E\/\k\. In fact the double change (k,E) — > (—k,—E) is 
equivalent to a time- inversion (see Figure I). 

Since x(t) is the quotient of two polynomials in t, with the denominator of higher degree than 
the numerator, the trajectories approach the origin when t increases. How can this behaviour be 
compatible with the conservation of the energy? An analysis of the expression that El takes on 
the trajectories shows that it reduces to the quotient of two functions both going down as t^oo, 
but in such a way that the ratio remains constant. Moreover it can also be proved that the velocity 
v x decreases as t^oo in such a way that the particle approaches but never reaches the origin in 
the phase plane. The important point is that the dependence of El with respect to v x is defined 
in such a way that even when dx/dt decreases the value of El remains constant. 

The phase space analysis shows that the origin is a nonelementary critical point for which the 
linear approximation is not valid. If we consider a small neighbourhood of the point, then we find, 
when k > 0, that it has four different sectors [22] (see Figure II): an elliptic sector (upper side), 
an hyperbolic sector (lower side) and two parabolic sectors (a source on the left and a sink on the 
right). 



3.2 Symplectic formalism and master symmetries 

It is known that the nonlinear equation (9) admits constants of motion depending explicitly on the 
time t [13], [16]. Now we prove that this property is a related with the fact that the Lagrangian (10) 
admits "master symmetries". This is in fact an important property, not only from the geometric 
point of view but also because it is directly related with the superintegrability of the two-dimensional 
version of this system. 

A function T that satisfies the following property 

d d m c jm+l 

dt ^ ' ' dt m r ' dt m + l 

is called a generator of integrals of motion of degree m. Notice that this means that the function 
T is a non-constant function generating a constant of motion by time derivation. If we denote by 
T x \ and T X 2, the functions 



T X 2 



v x + k x 2 ' v x + k x 2 

then we have 2 3 

Jt Tx2 = 1 > Jt Txl = kTx2 > d^ Txl = K d^ Txl = °- 
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Hence T x \ and T x2 are generators of integrals of motion for the Lagrangian (10). 

In geometric terms this property is related with the existence of "master symmetries" [23, 24, 
25, 26, 27]. If the dynamics is represented by a certain vector field T; then a vector field Z that 
satisfies the following two properties 

[z,r] = z/o, [z,r] = o, 

is called a "master symmetry" of degree m = 1 for V. If Z is such that 

[Z,T]=Z^0, [Z,T]^0 and [[Z,T],T}=0, 

then it is called a "master symmetry" of degree m = 2. 

We focus our attention in the case of master symmetries given rise via Z to constants of motion. 
Let L be a time-independent Lagrangian, Z the Hamiltonian vector field of a certain function T 
and suppose that Z is a (time- independent) master symmetry of T^. Then the time-dependent 
vector field Yz determined by Z and defined as 

Y z = Z + t[Z,T L ] + (^)t 2 [[Z,T L },T L } + ... 

is a time-dependent symmetry of I"l that, in the case m = 1, satisfies 

i(Yz)^ = d[T-tr L (T)] 

where = u)L+dEi/\dt. Hence the time-dependent function J* = T —tT l(T) is a time-dependent 
constant of motion (for m = 2 the corresponding constant J 1 is quadratic in t) 

We now return to the Lagrangian (10) and denote by Z x \ and Z x2 , the Hamiltonian vector fields 
of T xl and T x2 , 

i{Z x i)u) L = dT x i , i(Z x2 ) = dT x2 , 

which are given by 




where M x denotes M x = v x + k x 2 . Then we have 

[Z x i,T L ] = -kZ x2 , [[Z xl ,T L ] ,T L ] = 0. 

It is clear that Z x \ is a master symmetry. Concerning Z x2 it has some very interesting characteris- 
tics; it is a dynamical symmetry (that is, [Z x2 , T^] = 0) and it is a symplectic symmetry (that is, 
Cz x2 W£ = 0), but nevertheless it is not a Cartan symmetry of the Lagrangian system Tl because 
Z X 2(El)^0- Recall that the two properties Cx^L = and [X , Tl] = imply that CxdEL = 0, 
but from this we only obtain that X(El) must be a numerical constant. 
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As explained above, the vector fields Z x \ and Z X 2, determine two new vector fields Y\ and Y2, 
given by 

Y\ = Z X 2 + t [Z X 2, Tl] = Z X 2 , 

Y 2 = Z xl + t[Z xl ,T L ] + (^)t 2 [[Z xl ,T L ],T L ]=Z xl -ktZ X 2, 
such that they satisfy 

i{Yi)£l E = i(Z x2 )uJL + Z x2 {E L )dt = d[T x2 - t] , 

i(Y 2 ) n E = i(Z xl -kt Z x2 ) lo l + Z xl {E L ) dt - kt Z x2 {E L ) dt = d[T xl - ktT x2 + (^)kt 2 ] . 
Hence the two following functions 

J x i = T x2 -t, jl 2 =T xl -ktT X 2 + {^)kt 2 

are time-dependent integrals of motion determined by Z x \, Z X 2-, via Yi, Y2. 

In the next subsection we see that these symmetries and these time-dependent integrals are the 
origin of the n = 2 superintegrability. 



3.3 Two-dimensional nonlinear system: Lagrangian formalism and 
super-integrability 

We now want to study the nonlinear system 

representing the n = 2 version of the nonlinear equation (9). It is clear, from the results of n = 1, 
that these two equations can be considered as the Lagrange equations arising from the following 
two-dimensional Lagrangian 

L = \ , + . (12) 

v x + k\ x z v y + k 2 y z 

Therefore the dynamics is characterized by preserving the following Lagrangian energy 

(2v x + h x 2 ) (2v y + k 2 y 2 ) 



E L = - 



(v x + kix 2 ) 2 (v y + k2y 2 ) 2 



The first consequence of the Lagrangian character of the equations is that, as there is no coupling 
between the two degrees of freedom, the two one-dimensional energies are integrals of motion 

h _ (2^ + h x 2 ) ^ h _ (2v y + k2y 2 ) ^ SL Ii = Q i = 12 . 
(v x + ki x 2 ) 2 ' (v y + k 2 y 2 ) 2 dt 
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Our main goal in the study of this nonlinear problem is to prove that this system possesses 
the rather unusual property of super-integrability. At this point we recall that a system is called 
super-integrable if it is integrable in the Liouville-Arnold sense and, in addition, possesses more 
independent first integrals than degrees of freedom (see Refs. [28]- [36] for some articles published 
in these last years and Ref. [37] for a recent workshop on super-integrability). It is clear that, 
for this particular n = 2 system, super-integrability means the existence of a third independent 
integral I3 coupling the two degrees of freedom in similar way as the angular momentum for the 
isotropic linear harmonic oscillator. 



The four functions 



4x = T x2 -t, J t x2 = T xl -k l tT x2 + (^)k l t 2 , 
Jyi = T y2 -t, J t y2 =T yl -k 2 tT y2 + (±)k 2 t 2 , 



are time-dependent constants of motion 



— J*. - — 7* - i - 1 2 

dt dt yi ~ ' ' 

We can eliminate the time t by pairing these functions in two different ways and obtain the following 
two integrals of motion 

^3 = T x2 — T y2 , 

h = k 2 T x i + k{T y i - k\k 2 T x2 T y2 



that take the form 

/ 3 



x y 



v x + k\ x 2 v y + k 2 y 2 ' 



k 2 k\ h k 2 xy 

1 4 — — ; n + 



v x + kix 2 v y + k 2 y 2 {v x + kix 2 ){v y + k 2 y 2 )' 

Hence, the nonlinear system given by the equations (11) and characterizized by the Lagrangian 
(12) is a super integrable system. 

In geometric terms the symplectic form ul and the dynamical vector field are given by 

2dx/\dv x 2dyAdv y 
= 7 — ; — ^7 + 



(v x + kx 2 ) 3 (v y + k 2 y 2 ) 3 ' 
^ d d _, d _, d 

Tl = V X \-Vy — +F x h Fy — , 

ox oy ov x ovy 

where F x = — kix(3v x + k\x 2 ) and F y = —k 2 y(3v y + k 2 y 2 ). If we denote by Z xr and Z yr , the 
Hamiltonian vector fields of T xr and T yr , r = 1,2, then the vector field 

X3 = Z x2 — Zy 2 

is a dynamical symmetry, 

[x 3 ,r L ] = o, 
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as well a Cartan symmetry, 

X 3 (E L ) = 0, £x 3 ul = 0. 
It determines the function ^3 as the corresponding Hamiltonian 

i(X 3 )co L = dl 3 . 

Similarly the vector field X 4 given by 

X 4 = k 2 Z xl + k\Z yl - (k 1 k 2 ) (T y2 Z x2 + T x2 Z y2 ) 
is also a dynamical symmetry as well a Cartan symmetry, 

[X 4 , T L ] = , X 4 {E L ) = , C X4 oo L = 
and determines the function I 4 as the corresponding Hamiltonian 

i(X 4 ) 0JL = dl 4 . 

The main difference between these two symmetries is that X 3 is an exact symmetry of the 
Lagrangian, that is X 3 (L) = 0, and X 4 is a nonexact generalized Noether symmetry 

Ft = -fy£ + ± )t 

so that I4 is given by 

i{X 4 )6 L -F 4 = I 4 . 

An important property related with superintegrability is the existence of periodic orbits as, for 
example, in the Kepler problem or in the harmonic oscillator. This particular system is super- 
integrable but the motion, although bounded, is not periodic; instead of this the trajectories are 
"almost closed". We have found that the trajectories in the plane (x,y) are in fact "figure eight" 
curves with the particle making a complete circuit as t goes from — 00 to +00. The trajectory starts 
very close to the origin, passes through (0, 0) only once and returns once more to the origin, but 
for t — > —00 and t — > 00 only approaches (0, 0) as a limit. Figure III shows two curves in the plane 
(x,y) corresponding to two different values of E\ and E 2 . 



4 Lagrangian conservative approach to a nonlinear os- 
cillator 

4.1 One-dimensional system: Lagrangian formalism and integra- 
bility 

We now consider the following non-natural Lagrangian 

kv x + k 2 x 2 + w 2 ' ^ ^ 
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where k and w are arbitrary constants; the notation w 2 for the new parameter clearly advances 
that it will be interpreted as a frequency. We arrive to the following nonlinear equation 

(fix dx 

—r + 3kx{—) + k 2 x 3 + w 2 x = (14) 
at 1 at 

as well as to the following expression for the Lagrangian energy 

(2kv x + k 2 x 2 + w 2 ) 



E L = 



(kv x + k 2 x 2 + w 2 ) 2 



We can solve this new nonlinear problem by using the same technique as in the previous equation 
(9). Nevertheless we see that this new nonlinear equation is in fact a nonlinear oscillator and, 
because of this, we will use an approach that is as close as possible to the linear oscillator. 

We denote by X and W x the following two functions 

v x + kx 2 



X = -— = ^ , W, 



(kv x + k 2 x 2 + w 2 ) ' (k v x + k 2 x 2 + w 2 ) 



Then we have the following two properties: 
(i) The time evolution of X and W x is given by 

d „ „. d 



(ii) X and W x are related by 

w 2 X + k xW x = x . 

Next we see that property (i) is related with the conservation of the energy and property (ii) with 
the solution of the dynamics. In fact, from (i) we conclude that the function Ixw defined as 

Ixw = W 2 + w 2 X 2 . 

is an integral of the motion. We remark that, as the system is one-dimensional, the two integrals, 
El and Ixw, cannot be independent; in fact Ixw turns out to be the energy associate to an 
equivalent Lagrangian L of the form L = c\L + c^. A simple calculation gives 

E L = A(Z) - L = I xw with L = (^) 2 L-^. 



Also from (i) we arrive at 



^X + w 2 X = 0, ^ Wx + w 2 W x = 



so that X and W x are given by 



X = (— ) A sin (wt + (f>) , W x = A cos (wt + <f>) , A = \[E , 
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where E is the constant value of the new energy function El = Ixw- These two expressions, 
together with property (ii), lead to the following trigonometric function for the solution of the 
dynamics 

w \fE sin (wt + <fi) 

x = = . 

1 - k\/E cos (wt + 4>) 

Figure IV represents x(t) as a function of t for several values of the energy E in the regular 
oscillatory region (0 < E < 1/k 2 ). It is clear that for small values of E the oscillations are rather 
similar to the oscillations of the linear system, but for other values of E (roughly speaking, for 
E > 0.3/fc 2 ) the nonlinearity introduces deformations in the oscillations. For higher values of E 
the nonlinearity changes drastically the aspect of the solution. Figure V shows the phase portrait; 
it clearly shows that the motions from left to right and from right to left take place at different 
velocities. 



4.2 Two-dimensional system: Superintegrability and nonlinear 
Lissajous figures 

We now turn our attention to the two-dimensional version of this nonlinear oscillator. If we consider 
the following Lagrangian 

1 1 

^ k\ v x + k\ x 2 + w\ ~*~ k 2 v y + k 2 y 2 + w\ 

then we arrive at the following equations 

d 2 x „ , , dx 



+ 3k 1 x(—) + k(x 6 + w(x = 0, 



dt 2 K dt 

+3k 2 yA + k 2 y* + w 2 y = (16) 



dt 2 zyy dt 
and to the following expression for the energy function 

E _ _ (2fci Vx + k\x 2 + wf) _ (2 k 2 v y + k\y 2 + w\) 
(k\ v x + k\ x 2 + w\ ) 2 (v y + k 2 y 2 + w^) 2 

In this n = 2 case we have two pairs of functions (Xj,Wj), j = 1, 2, 

v X rxr Vx + klX 2 

A l — —P) 9 5 W l 



k\ v x + k\x 2 + w\ ' k\ v x + k\ x 2 + w 2 ' 



in such a way that we have 



k 2 v y + k\y 2 + W 2 ,' k 2 v y + k%y 2 + w\ 



4 W j = - w 2 Xj , ^-Xj = Wj 
dt 1 3 J ' dt 3 3 
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and 2 

^2 + ™f ^ = > ^2 ^ + ™f ^ = • 

A similar calculation shows that 

w\ X 1 w\ X 2 

x — .. — - 



l-JfeiW^i' 1 - fc 2 W 2 ' 

from which we obtain solution of the dynamics 

w\ y/Ei sin (w\t + (pi) w 2 \fE~2~ sin (u;2i + 4>2) 



x 



1 — ki \[E[ cos (u;it + (pi) ' 1 — \fE~2 cos (u^i + ^2) 



We now study the superintegrability of the rational case, that is, wi = niwo and w 2 = n 2 wo, 
with m and n 2 positive integral numbers. 

Proposition 2 Let Ki and K2 be the following two functions 

Kj = Wj+i njWo Xj , j = 1,2. 
XTien the complex functions Kjj defined as 

Hi, = (K*) n * , i,j = 1,2, 

are constants of the motion. 

Proof: The time evolution of the functions K4 and K 2 is given by 

d d d 

On the other side we have 

and from here the property follows by direct calculus. 
Thus the three functions 

/ 1 = |Ki| 2 , / 2 = |K 2 | 2 , Ii 2 = K™(K*) n \ 

are constants of the motion. The two first functions, Ii and I 2 , are the two one-degree of freedom 
energies; concerning Ii 2 , as it has complex value 

hi = h + \h, 

it determines not just one but two real constants of motion. Of course, if we consider ^3 as the 
new additional constant, I4 is a function of ii, I 2 and I3. We thus conclude that the existence of 
superintegrabilty and periodic trajectories (Lissajous figures) is preserved by the nonlinearity. 
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We particularize these results for the two first commensurable cases. In the isotropic case, 
wi = w 2 = wq, the two functions, ^3 and I4, are given by 



h = x 1 w 2 -x 2 w l 



(xvy - yv x ) + (k 2 y - k x x)xy 
(ki v x + k\ x 2 + wl)(k 2 v y + k 2 y 2 + Wq) ' 



h = w lW2 + w \x x x 2 - {v * + kl x2){v y ± h y2) + w ° xy 



(fei + fef X 2 + WJq)(^2 % + ^2 V 2 + W o) ' 

representing the nonlinear versions of the angular momentum and the nondiagonal component of 
the Fradkin tensor respectively [38, 39]. In fact, when k±,k 2 — > 0, these two function reduce to the 
appropriate expressions 

lim fc ^o^3 = ( \ ) (xv y - yv x ) , 

lim ^0-^4 = (— ){v x v y + wlxy) . 
w Q 

Of course, for wq — ► 0, we recover the ^3 obtained in the previous section; I4 just reduces to a 
trivial numerical constant. 



lim w ^oh 
lim w ^oU 



{xvy - yv x ) + {k 2 y - k\_x)xy 

k\k 2 {v x + ki x 2 )(v y + k 2 y 2 ) 
1 



k\k 2 

Figure VI represents some closed trajectories in the plane (x,y); it is clear that for small energies 
the curves look rather similar to the ellipses of the linear case, but for other (not so small) values 
of E the curves lose their elliptic shape and adopt other not so symmetric forms. 

Now, we consider the anisotropic case w\ = wq, w 2 = 2wq. Then I 3 and I4 are given by 

h = {x x w 2 -x 2 w x )w x + wlxlx 2 

(v x + hx 2 ) [ (xv y - yv x ) + (k 2 y - hx)xy } + w^x 2 y 
(k lVx + k 2 x 2 + wl) 2 {k 2 v y + kly 2 + 4wg) 
I 4 = wfW 2 + wl{AX 2 Wi - X l W 2 )X 1 

(v x + ki x 2 ) 2 (v y + k 2 y 2 ) + wl [4yv x - xv y + (4^x - k 2 y)xy]x 



{ki v x + k\ x 2 + wl) 2 (k 2 v y + k^y 2 + Aw^) 

Figure VII represents two nonlinear Lissajous figures in the plane (x,y). The situation is similar 
to that of figure VI; close resemblance with the linear figures for small values of the energies and 
rather strange figures for other values of E. We must say that, in this case, the form is strongly 
dependent on the phase difference <t>i 2 = 4>i — 4> 2 - 



5 Final Comments 

We have studied two nonlinear systems using, as starting point, the important property of the 
Lagrangian origin of the second-order Riccati equations. In this way we could use the constant value 
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E of the energy as an appropriate parameter for characterizing the behaviour of these two systems. 
Moreover we have proved the existence of two-dimensional versions endowed with superintegrability 
and we have obtained the explicit expressions of the additional integrals. 

Concerning the superintegrability, we recall that most of known superintegrable systems are su- 
perseparable systems, that is, systems that admit Hamilton-Jacobi separation of variables (Schro- 
dinger in the quantum case) in more than one coordinate system. Nevertheless as all the systems 
studied in this paper are nonlinear systems with a nonstandard Lagrangian the superintegrabil- 
ity has been proved by considering a different approach and without making use of the multiple 
separability. In spite of this the possible separation of variables in other coordinate systems than 
Cartesian ones must be studied. 

We also mention that these two nonlinear systems can be generalized in several different ways 
by making use of the Lagrangian approach. Firstly, we recall that we have proved in Sec. 2 the 
Lagrangian origin not only of the Riccati systems but also of the more general equation (4); this 
means that these more general equations could also be studied by making use of the conservation 
of the energy. Secondly, the Lagrangian (6) admits the following natural generalization 

! (17) 



m(x)v x + kU (x) 



that leads to a nonlinear second-order equation with a position-dependent effective mass and a ad- 
ditional dissipative-looking term of the form m'(x)v 2 . This new Lagrangian (17) seems interesting, 
not only because it generalizes (6) but also because it has a more direct geometrical interpretation 
since the linear function m(x)v x can be considered as associated to the one-form \i = m(x)dx. 

The Lagrangian (6) has, as associated Hamiltonian, the following function 

H = -2^=p~ x -kU{x)p x (18) 

that, in addition to its nonnatural character (as was to be expected), has the annoying presence of 
the momentum inside a root. In the particular case of the nonlinear oscillator the Hamiltonian is 
given as follows 

■2W , : „ o W 2 



H = — 



V-kp x + (kx 2 + ^-) Px | + ^ . (19) 



Nevertheless in this case we have an important property; in fact, if we make use of the canonical 
transformation (x,p x ) — > (Q,P) given by 

y/2 \[2 
Q = (—) x V-kpx , P = (—) I 1 ~ wy/-kp x ] , 



H=(l)(P 2 + w 2 Q 2 ) 



then we arrive at 

H=( 

Hence the very peculiar Hamiltonian (19) and the standard linear oscillator are canonically related. 
Nevertheless note that this transformation has a nonpoint character and, because of this, it cannot 
be directly used in the Lagrangian approach. This brings up the question of the possible existence 
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of similar nonpoint transformations for Hamiltonians obtained from other Riccati Lagrangians. We 
think that this possibility is a open question to be studied. 

Finally, we mention the study of the quantized versions of all these nonlinear systems. We note 
that this question must be carry out only after the obtaining of the appropriate Hamiltonian versions 
(the direct quantum study of the Lagrangian equations appears as a difficult task). Nevertheless 
as the Lagrangians are nonstandard the Hamiltonians also appear with a unusual dependence of 
the momenta (see the above expressions). In any case the possibility or impossibility of quantizing 
these systems is a matter that must be investigated. 
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ure Captions 



Figure I. Plot of x as a function of t, for k = 1 and three different values of the energy: 
E = -0.5, E = -1.0, and E = -1.5. 

Figure II. Phase space trajectories, for fc = 1, in the neighbourhood of the origin. 

Figure III. "Figure eight" in the plane (x, y) corresponding to k± = &2 = 1 and energies 
Ei = — 1, E2 = — 5 (thick curve) and E\ = — 1, E2 = —10 (dash curve). 

Figure IV. Plot of x as a function of t, for (A; = 1, w = 1) and three different values of the 
energy: E = 0.2 (small thick curve), E = 0.5 (middle curve), and E = 0.8 (curve with great 
oscillations). The E = 0.2 curve is very similar to a pure sine or cosine curve but, for higher 
values of E, the plot shows clearly the effects of the nonlinearity. 

Figure V. Phase trajectories corresponding to four different values of the energy (E = 
0.2,0.4,0.6, and 0.8). The trajectories are closed curves representing periodic motions that 
for small values of E can be considered as rather similar to ellipses; for other values of E the 
curves modify their shape and they lengthen towards the upper side of the phase plane. The 
motion is asymmetric in the sense that the particle moves from right to left in a slowly way 
but returns, from left to right, with a much higher velocity that takes its maximum value at 
the center point x = 0. 

Figure VI. Closed trajectories (ellipses in the linear case) in the plane (x, y) corresponding 
to W\ = 1, W2 = 1 and four different values of the energy (E\ = £2 = 0.2,0.4,0.6, and 0.8). 

Figure VII. Nonlinear Lissajous figures in the plane (x, y) corresponding to "figure eight" 
trajectories associated to W\ = 1, wi = 2 and energies E\ = E2 = 0.2 (small eight-looking 
curve) and E\ = E2 = 0.6 (big butterfly-looking curve). 
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Figure I. Plot of x as a function of t, for k = 1 and three different values of the energy: 
E = -0.5, E = -1.0, and E = -1.5. 



Figure II. Phase space trajectories, for k = 1, in the neighbourhood of the origin. 
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Figure III. "Figure eight" in the plane (x, y) corresponding to k\ = k 2 = 1 and energies 
E\ = — 1,E 2 = —5 (thick curve) and E\ = — 1,E 2 = —10 (dash curve). 




Figure IV. Plot of x as a function of t, for (k = l,w = 1) and three different values of the 
energy: E = 0.2 (small thick curve), E = 0.5 (middle curve), and E = 0.8 (curve with great 
oscillations). The E = 0.2 curve is very similar to a pure sine or cosine curve but, for higher values 
of E, the plot shows clearly the effects of the nonlinearity. 
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Figure V. Phase trajectories corresponding to four different values of the energy (E = 0.2, 0.4, 0.6, 
and 0.8). The trajectories are closed curves representing periodic motions that for small values of 
E can be considered as rather similar to ellipses; for other values of E the curves modify their 
shape and they lengthen towards the upper side of the phase plane. The motion is asymmetric in 
the sense that the particle moves from right to left in a slowly way but returns, from left to right, 
with a much higher velocity that takes its maximum value at the center point x = 0. 
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Figure VII. Nonlinear Lissajous figures in the plane (x,y) corresponding to "figure eight" 
trajectories associated to w\ = 1, wi = 2 and energies E\ = E2 = 0.2 (small eight-looking curve) 
and Ei = E2 = 0.6 (big butterfly-looking curve). 
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